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ABSTRACT 



The results of an experimental investigation of the effect 
of adverse pressure gradient on the vortex breakdown phenom- 
enon in a diverging tube are presented. Adverse pressure 
gradient was found to be as significant in determining the 
breakdown position as were the previously known parameters, 
namely, the Reynolds and circulation numbers. It was found 
that an increase in any of the three parameters serves to 
move the breakdown upstream, toward the origin of the vortex. 
It was further found that beyond some range, an increase in 
adverse pressure gradient causes boundary-layer separation 
in the diverging tube. Separation restricts the swirling 
flow and thus limits the effective adverse pressure gradient 
acting on the vortex. Data are presented to illustrate 
these effects. An approximate momentum analysis for pre- 
dicting the vortex breakdown position was carried out. The 
results were found to be only marginally satisfactory. In 
spite of that, however, the method appears to hold much 
promise as a breakdown predictor. 
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I . INTRODUCTION 



Vortex breakdown is an abrupt change in the structure 
of the core of a swirling flow. It was first noted by Peck- 
ham and Atkinson in 1957 during their studies of air flow 
over delta wings at large angles of attack [1]. Since then 
vortex breakdowns have also been observed in straight or 
diverging cylindrical tubes, cyclone separators, behind tor- 
pedos, and in both nozzels and diffusers. Vortex breakdown 
may in fact occur in any vortex if the critical swirl and 
flow conditions are met. Hence the study of this phenomenon 
may have a practical application in any device containing a 
swirling fluid. 

Three distinct types of vortex breakdown have been re- 
ported by Sarpkaya [2], These are the spiral, axisymmetric 
and double-helix types. Both the spiral and axisymmetric 
types initially form at a stagnation point along the vortex 
core, or centerline. These breakdowns mark the transition 
between laminar and turbulent flows. 

In the case of spiral breakdown, the core appears to 
first suddenly stagnate, then kink to one side, and then 
form a spiral as the kink rotates. After a few revolutions 
the spiral becomes lost in the violently turbulent fluid 
flow. 

The axisymmetric breakdown also initially appears as 
a stagnation in the flow, followed by a roughly right-angled 
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kink. However, rather than forming a spiral, the rotating 
kink forms the forward boundary of an egg-shaped body, which 
fills from the rear, exhibits significant internal backflow 
and turbulence, and then discharges through the downstream 
side of the body. Depending on the magnitudes of both the 
circulation and flow rate (i.e., Reynolds number) in the 
vortex, the axisymmetric body may be closed from as little 
as 50% of its projected length, 1 up to 95%. When the body 
is only partially closed, turbulence commences at the exit 
of the body, and fans out to the outer reaches of the still 
laminar vortex in a very rough conical shape. When the 
axisymmetric "bubble" is almost fully closed, the exit core 
appears to re-form for one or two body lengths, and then 
either kink again and form a spiral as described above or 
proceed immediately to full turbulence in a smooth conical 
shape . 

The double-helix type of breakdown does not include a 
definite stagnation point in the core; rather, a mild dis- 
turbance causes a slight slowing and swelling of the core. 
The core then shears into a tape-like helix formation. Its 
appearance is that of two ribbons emanating from a slight 
expansion of the core streakline, then entwining about each 



1 This type is sometimes referred to as a "near- 
axismmetric" type [2]. 
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other. In the case of the double-helix breakdown, which 
generally occur at low Reynolds numbers (Re ~ 700 to 1500), 
there is no immediate transition from laminar to turbulent 
flow. It may therefore be considered by some to be more of 
a flow instability, significantly different in nature from 
the more common and relatively abrupt spiral and axisymmetric 
vortex breakdown forms [2]. 

Shortly after Peckham and Atkinson reported their dis- 
covery, numerous researchers attempted to develop a suitable 
model to describe and predict the occurence of vortex break- 
down. The experimental contributions of the various workers 
have been well documented by Robertson [3], and will not be 
pursued herein. Alongside the experimental efforts, several 
theoretical analyses were developed, but the divergence of 
approach and initial assumptions have created significantly 
different and conflicting theories to explain the vortex 
breakdown phenomenon. Recently Hall [4] presented an excel- 
lent summary of these studies and compared the various anal- 
yses, highlighting their major strengths and weaknesses. 

He grouped the various proposals into three categories, in 
which the basic ideas were as follows: 

1. The phenomenon is in some sense like the separation 
of a two-dimensional boundary layer (Gart shore 
1962, Hall 1967 , and others). 

2. The phenomenon is a consequence of hydrodynamic 
instability (Ludwieg 1962, 1965 , and others). 

3. The phenomenon depends in an essential way on 
the existance of a critical state (Squire I960, 
Benjamin 1962, 1967, Bossel 1967> 1969)* The 
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proposals in this category have equivalent 
definitions for the critical state but a 
range of interpretations of its significance. 

However, after detailed analysis of the theories. Hall con- 
cluded that none of the proposals was adequate to completely 
explain what a vortex breakdown was, and when or where it 
would occur. 

Interestingly enough, Sarpkaya [5] noted that while the 
axisymmetric vortex breakdown basically corresponded to the 
finite-transition critical-state theory of Benjamin [6], 
the double-helix and spiral types were more in accord with 
Ludwieg's instability theory [73* Carrying this one step 
further, Mager [8,9] presented a mathematical analysis based 
on the works of Hall and Sarpkaya, which essentially com- 
bined the finite-transition and hydrodynamic instability 
theories. In essence he showed that each theory was basi- 
cally correct as far as it went, but that each was only a 
portion of a larger theory. He proposed tha't a complete 
analysis should, not stop when the breakdown occurred, but 
rather carry through from the laminar vortex to the onset of 
turbulent mixing [9]. This proposal will be discussed further. 

Prom the foregoing, it should be obvious that much is 
still not known about the vortex breakdown phenomenon. Not 
only is there controversy over how or why the breakdown 
occurs, but also on the degree of importance of the various 
parameters; e.g., velocity profile, magnitude of swirl, Rey- 
nolds number and the effect of axial pressure gradients. 
Sarpkaya [2] has shown that the first three of these are 
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quite significant, at least in the study of breakdowns in 
water. Several researchers including Hall [4], Jones [10], 
Lambourne and Bryer [11], Mager [8] and Sarpkaya [2], have 
concluded that an adverse (positive) axial pressure gradient 
in the flow has an effect on the occurrence of breakdowns. 
Lowson [12], as a result of his analysis and experiments, 
reported that "the final position of the majority of the 
practical breakdowns is due primarily to the reaction of 
pressure gradient effects on an already established break- 
down." However, as yet, no one has definitely verified either 
theoretically or experimentally the role that the axial 
pressure gradient is presumed to play. 

The purpose of the present investigation then was two- 
fold. The first objective was to examine what, if any, 
effect adverse pressure gradient has' on vortex breakdowns. 

This was done by using cylindrical test sections of varying 
degrees of divergence to provide different pressure gradients 
to the vortex flow. The second objective was to briefly 
examine Mager* s theoretical analysis of vortex breakdown, 
and to determine whether it offered any more promise as a 
vortex breakdown predictor than its predecessors. Considerable 
data were taken on vortex breakdown position and type over 
a wide range of flow parameters, and this was used to both 
graphically demonstrate the effect of the adverse pressure 
gradient and provide data for a computer analysis of Mager *s 
hypothesis as suitably adapted to the present experimental 
system. 



13 



II. DESCRIPTION OF THE EXPERIMENTAL APPARATUS 



The experimental apparatus consisted of a test tank, 
water supply reservoir, centrifugal pump, overhead reservoir, 
inlet box, two dye injection lines, vane-angle adjustment 
mechanism, thermometer, vortex-breakdown posit ion/swirl 
angle measuring device, two rotameters and necessary piping 
and valves, connected as shown in Figure 1. 

A scale drawing of the test tank is presented in Figure 
2. It was constructed of 3/^" Plexiglas plates, held to- 
gether by screws and sealed with silicon sealant. The re- 
movable top was made in two parts for ease of handling during 
removal for changing the test sections, and was fitted to 
the body with a silicon rubber gasket. Internal fittings in 
the test tank were as follows: 



1. Sectioned baffle-plate - This was made- of 3/^" 
Plexiglas, and drilled with 39 rows of holes 
stepped from 1/8" at the bottom to 5/16" on 
top. This provided a nearly uniform velocity 
distribution . in the test tank. Fine mesh 
screens were installed on each side of the 
baffle-plate sections to reduce the turbulence 
of the entering fluid. 

2 . Adjustable vane mechanism - An additional 3/^" 
Plexiglas plate was attached to the upstream 
wall of the test tank. A movable circular ring, 
9" in diameter, was set into this plate. 
Thirty-two fixed pins were attached to the ring 
in a symmetric circular array, 11.25 degrees 
apart. Each of these pins fitted into a slot 

in a streamlined blade, or vane. The vanes 
were pinned to the inner plate, and pivoted 
together about their fixed pins as the movable 
ring was rotated by a hand-screw mechanism. 

The indicator on the hand-screw mechanism was 



calibrated to read the vane angle directly to 
the nearest half degree, and was capable of 
changing vane angle from plus to minus sixty 
degrees, thus varying the amount of swirl 
imparted to the vortex. Blade dimensions and 
geometry are shown in Figure 3. 

3. Cone and bell-mouth - A six inch diameter cone 
was attached to the additional plate described 
above, inside the circle of adjustable blades. 

This cone and its matching bell-mouth formed 
the entrance to the test section. The center- 
line was bored to provide an entrance for the core 
dye stream. An outer dye stream (to show the swirl 
angle) was introduced through an adjustable probe 
located in the flow stream between the core and 
bell-mouth. To ensure that the bell-mouth (and 
the connecting test section) were centered o s n 

the cone, simple hand-screw adjustments were 
provided for both horizontal and vertical 
centering. 

4 . Divergent test section and outlet tube - Four 
pairs of these were made in matched sets, machined 
to close tolerances to ensure that there would 

be no obstruction to smooth flow through the bell- 
mouth, test section and outlet tube. To test 
the effect of pressure gradient, tubes of varying 
angles of divergence were manufactured and in- 
stalled. The inlet diameter of all test sections 
was identical to that of the bell-mouth outlet, 
1.755".’ Test section # 1 (TS-1) had a length 
of 10 3/8", and an outlet diameter of 2 1/4", 
providing a (half) angle of divergence (y) of 
1.38°. Test section #2 was made by adding a five 
inch extension to TS-1, maintaining the same angle 
of divergence, thereby giving it an outlet di- 
ameter of 2 1/2". TS-3 had a length of 15", and 
y = 2 . 36 °; and, TS-4 had a length of 15 3/8", 
and y = 3-25°. Each test section was installed 
with an outlet tube that matched the test sec- 
tion diameter. The outlet tube was mated with 
a matching nozzle at its far end, to constrict 
the flow into the 1 1/2" outlet pipe leading 
to the rotameter bank. 



Exterior fittings on the test tank were the inlet box, 
thermometer to measure water temperature, and a 26" over- 
flow column to maintain constant tank pressure and absorb 
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flow fluctuations. A small device to measure breakdown 
position and swirl angle rolled along a section of model 
railroad track permanently attached to the top of the test 
tank. A one inch deep Plexiglas sighting bar was suspen- 
ded from a wheeled, counter-weighted cart. The sighting 
bar also provided a suitable location to mount a fixed pro- 
tractor and rotatable straightedge to measure swirl angle. 

To keep the system as free of disturbances as possible, 
air vents were provided at several locations, and the entire 
test tank was mounted on foam rubber. 

Plow of water through the system can readily be traced 
by referring to Figure 1. City water entered the lower 
reservoir, was pumped to the overhead reservoir, flowed by 
gravity into the inlet box, and entered the test tank. It 
then passed through the baffle-plate and screens, slowly 
traveled the length of the test tank, passed through the 
adjustable vanes, and between the cone and be II -mouth where 
it formed a laminar vortex. The water then entered the 
diverging test section and exited the tank through the outlet 
tube. After passing through the outlet nozzle, the water 
was carried to either (or both) of the rotameters, and 
discharged to a drain. Rotameter ratings at 100% flow were 
10.0 and 2.70 gallons per minute respectively, and provided 
the capability of achieving Reynolds numbers within a range 
of Re = 300 to 18,000. 

A photograph of the system is presented in Figure 4. 

Out of sight are the pump, lower reservoir, inlet box and 
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the overhead reservoir. The large cement column to the 
rear of the system supports the overhead reservoir, located 
approximately fifteen feet above the test tank. 
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III. EXPERIMENTAL PROCEDURE 



The initial phase of this investigation consisted of 
determining the location of vortex breakdowns in tubes of 
various angles of divergence for identical flow conditions. 
The variables which could be controlled were circulation, 
Reynolds number, and degree of divergence of the test sec- 
tion. Circulation (ft), a measure of the vortex strength, 
was taken as the product of the tangential component of the 
velocity of the fluid leaving the vanes and the circumfer- 
ence of a circle drawn through the vane tips. Prom conti- 
nuity and geometry, this resulted in the following equation 
C2], 



f = 2ir(R v - m cos 4>) V tQtal sin \p 



where <f> is the vane angle (read directly from the indicator 
in degrees) and the other terms are calculated or defined 
as shown in Figure 3. Circulation was normalized by defining 
a circulation number as 



ft = r/u fc D fc 



which may be expressed as 



ft = 



tt D. (R - m cos <j)) sin 

X> V 

64 d h 
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A simple computer program was written to calculate the 
corresponding circulation number for each vane angle 
setting. 

Reynolds number (Re) was used to normalize the flow 
rate through the test section. 



where was the test section diameter at the start of 
divergence (1.755"), U^. was uniform velocity, and v the 
kinematic viscosity of the working fluid, water. was 

determined from the flow rate measured with a calibrated 
rotameter. As the experiment was run using city water, v 
did not remain constant, but varied as a function of water 
temperature. This variation has been taken into considera- 
tion in the calculation of Reynolds number by using the 
appropriate value of kinematic viscosity for the temperature 
measured. 

The experimental effort was conducted in four parts 
(one for each test section), and each part was conducted in 
an identical manner. The test tank was filled very slowly 
from the overhead reservoir in order to minimize turbulence 
in the test tank. Air vents, located at strategic points, 
were opened to eliminate any air pockets that might other- 
wise have caused disturbances in the breakdown position. 
After the test tank was filled, it was allowed to remain 
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idle for approximately one hour to permit the remaining 

2 

small amount of turbulence to dissipate. When the fluid 
was presumed still, the water supply was resumed and water 
run at a low flow rate through the entire system. 

Data runs were made by selecting a suitable vane angle 
and holding it constant while varying the flow rate through 
the test section by carefully setting the rotameter valve 
or valves. Vane angle setting varied from every two or 
three degrees for lower circulation numbers to five degrees 
for higher circulations. The total range of circulation 
numbers for all test sections varied from ft = 0.50 to 5*00 
(<{> = 13° to 55° ), although only a portion of this range 
could be used for each test section. 

For a fixed vane angle, the rotameter setting was in- 
creased in small increments, and conditions in the test 
section recorded. When a vortex breakdown was present, the 
following information was recorded: 

1. Vane angle (° - to nearest 0.2°) 

2. Flow rate {% - to nearest 0.5$) 

3. Water temperature (°F - to nearest 0.2°) 

■4. Swirl angle^ (° - to nearest 2°) 

5. Breakdown type, or types if in hysteresis region 

6. Breakdown position -(inches from start of test 
section divergence - to nearest 0.1") 



It would have been preferable to allow the water to 
remain in the test tank from day to day, but this would have 
permitted dissolved air to settle throughout the system and 
disturb the flow. 

J Swirl angle was taken as the angle of the tangent to 
the outer dye stream at the nearest crossing of the core 
prior to the breakdown. It is considered by some to be an 
indication of a critical flow condition in the vortex [2], 
but it was not used in this investigation. 
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7. Breakdown diameter (inches - to nearest 1/16") 

8. Percent that breakdown was closed, if axisym- 
metric type ($ - to nearest 5$) 

9. Ratio of length to diameter, if axisymmetric 
type ($ - to nearest 25$) 

10. Comments - anything peculiar to the current 
data point 

Breakdown and swirl visualization was accomplished by 
injecting dye into the vortex as previously described. 
Commercial food coloring at full strength was used as the 
dye , and although it was known to be slightly denser than 
water, the effect of the density difference was not 
noticeable . 

The most time consuming portion of the experiment was 
marking the location of the vortex breakdown. As the break- 
down is a dynamic phenomenon, its position oscillated within 
a narrow range, even though circulation and flow rate re- 
mained constant. All data points taken exhibited a slight 
random fluctuation in the axial direction. The average 
distance traveled during these fluctuations was approximately 
± 0.2 in. As experience was gained, it was found that by 
following the slow motion of the breakdown with the position 
indicator for about three minutes and mentally noting the 
limits, there would be roughly a 95$ probability that the 
breakdown would remain in a 0.4" range. (This was true for 
both large and small flow rates.) The midpoint of this 
range was taken to be the (average) breakdown position. 

This procedure was further complicated in the hysteresis 
region when two (or more) breakdown types were alternately 
present for the identical flow and circulation conditions. 
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When this occurred, it was necessary not only to determine 
the average position for each, but also the relative propor- 
tion of time that each was present. This latter step was 
necessary to facilitate the preparation of graphs in which 
average breakdown position for a given flow and circulation 
setting would be shown. 

Finally, after two complete sets of data runs were made, 
photographs were taken of various breakdowns over a wide 
range of flow conditions. These were taken on 2 1/4" x 
2 1/4" negatives in both color and black and white, using 
only the lights installed to the rear of, and under the 
test tank. Several of these photographs have been reproduced 
as Figures 5 through 13, and 19 through 23. 

Upon completion of data runs and photography, the test 
tank was opened and cleaned, and a different test section 
fitted and installed. As the test sections were partially 
completed in advance, this step was accomplished within a 
day or two. 
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IV. ANALYTICAL CONSIDERATIONS 

An exploratory approximate analysis of the effect of 
adverse pressure gradient on the location of vortex break- 
down consisted of the consideration of the equations of 
conservation of mass and momentum, the representation of 
the approximate velocity profiles by suitable polynomials, 
and the calculation of the distance at which the axial 
velocity gradient increases rapidly leading to a stagnation 
point . 

The equations of motion for axially symmetric, incom- 
pressible flow, with gradients in the axial direction much 
smaller than those in the radial direction, may be written 
as [8] 



(ur) + (wr) z 




0 



( 1 ) 

( 2 ) 



(uvr 2 ) r + (vwr 2 ) z = j^-[r 2 (v r - J)] r (3) 



(uwr) + 
r 



, 2 \ 
(w r) z 



_(H) P + -L 

K 2 )r z. + Re 



a 



•[r (w) ] 

u \ ' p J p 



( 4 ) 



where Re = Q R./v, and all velocities are normalized by 

3 . X/ 

the maximum velocity Q, all pressures normalized by the 

_2 

maximum dynamic pressure 0.5pQ , and all lengths normalized 
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by the radius of the tube at the start of the diverging 
test section . 

The above equations can be reduced to the following 
integral equations [8] 



where 



and 



with 



X 12 ' = 2/Re (5) 

I,-,’ + % I ?? ’ + I,W» - (S/2) 2 *(a’/a) = 0 

11 22 1 ( 6 ) 



6 

1^ E / (W - w)r dr, 
o 

1-^ = / w(W - w)r dr, 

o 

6 

I 12 E / w[ 1 - (^)]r dr, 
o 

6 2 

I 22 E / v r dr, 
o 



(7) 



P ’ = 2(w w ' - 
o o o 



S % 2 



WW’ ) + (-) 
a 



- (Ap)' = 



Rk (W rr } o 

( 8 ) 




6 2 
V 

/ — dr 
r 



(9) 



2 _ 

where a = 6 and S = T /R, Q are the normalized core area 

a t 

and the normalized circulation respectively. 



The foregoing equations may be integrated further only 

if one assumes a set of suitable velocity profiles. This 

procedure is quite similar to that commonly used in employing 

4 

the integral-momentum equation of Karman. One may do so 
by setting 



w - W[a + (1 - a) f^, v = V(f 2 + 3^) 



V = r a /6 



( 10 ) 



where 



n = r/6 , 

f 2 = n(2 - n 2 ) 



f x = n 2 (6 - 8n + 3n 2 ) 
f ^ = nd - n ) 2 



(ID 



These velocity profiles satisfy the boundary conditions at 
r = 0 and r = <$ , and reduce to the swirling-potential-flow 
velocity distribution outside the vortex core. The parameters 
a and 3 are initially unknown and will be determined later. 
Finally, to complete the formulation of the analysis, one 
needs to express the conservation of massflow within the 
diverging test section by writing 

M = WA - 21 _ (12) 

^For an example of this procedure, see Schlicting [133- 
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where A is given by 



A = (1 + Z tan y) 2 . (13) 

Z is the axial distance, and y is the half-angle of test 
section divergence in Equation (13). Evaluating the inte- 
grals given by Equation (7) by the use of the velocity pro- 
files, Equations (10) and (11), and inserting the various 
I. , values in Equations (5, 6, 8, 9, and 12), produces the 
following four equations: 

[w 2 (K 10 + aK 1]L + K 12 a 2 ) - r 2 /4a]a' +[aW 2 (K i:L + 2aK 12 )]a' 
+ [%r 2 (b 1]L + 2Bb 12 )]3 l + [2aW(K 1() + K^a + K 12 ct 2 ) 

+ 0.1 aW(l - a)]W' = [0] (14) 

[-0 . 2W(1 - a)]a' + [0 . 2aW]a ' + [A - 0.2a(l - cc)]W' 

= [-WA* ] (15) 

[ (T/a) 2 (1 + b 20 + b 21 3 + b 22 3 2 )]a' + [2ctW 2 ]a' 

+ [-(b 21 + 23b 22 ) r 2 /a]e' + [2a 2 W - 2W]W' 

= [48(1 - a )/(a Re)] (16) 

[W(K 20 + aK 21 )]a' + [aWK 21 ]a' + [aW(b 31 + ab 1|1 )]3’ 

+ [a(K 2Q + aK 21 )] = [2/Re] (17) 
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